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ABSTRACT 

Since cosmic-microwave-background (CMB) and large-scale-structure (LSS) data will 
shortly improve dramatically with the Microwave Anisotropy Probe (MAP) and Planck Sur- 
veyor, and the Anglo-Australian 2-Degree Field (2dF) and Sloan Digital Sky Survey (SDSS), 
respectively, it is timely to ask which of the CMB or LSS will provide a better probe of pri- 
mordial non-gaussianity. In this paper we consider this question, using the bispectrum as a 
discriminating statistic. We consider several non-gaussian models and find that in each case 
the CMB will provide a better probe of primordial non-gaussianity. Since the bispectrum is the 
lowest-order statistic expected to arise in a generic non-gaussian model, our results suggest 
that if CMB maps appear gaussian, then apparent deviations from gaussian initial conditions 
in galaxy surveys can be attributed with confidence to the effects of biasing. We demonstrate 
this precisely for the spatial bispectrum induced by local non-linear biasing. 



1 INTRODUCTION 



It is widely accepted that the large-scale structures we observe in 
the Universe today originated from gravitational evolution of small 
primordial fluctuations in the matter density. Information about the 
physical processes that generated these primordial fluctuations can 
be gleaned by testing whether their statistical distribution is well 
approximated by a gaussian random field. In particular, the sim- 
plest versions of inflation predict gaussian initial fluctuations (e.g., 
Guth & Pi 1982; Hawking 1982; Starobinsky 1982; Bardeen, Stein- 
hardt & Turner 1983), but there are other models of inflation (Allen, 
Grinstein & Wise 1987; Kofman & Pogosian 1988; Salopek, Bond 
& Bardeen 1989) and models where structure is seeded by topo- 
logical defects (Vilenkin 1985; Vachaspati 1986; Hill, Schramm 
& Fry 1989; Turok 1989; Albrecht & Stebbins 1992) that gener- 
ate non-gaussian fluctuations. By looking at cosmic-microwave- 
background (CMB) anisotropies we can probe cosmic fluctuations 
at a time when their statistical distribution should have been close 
to their original form. At present, the limited signal-to-noise or sky 
coverage of existing experiments is not sufficient to provide con- 
clusive evidence either for or against non-gaussianity (e.g., Heav- 
ens 1998; Ferreira, Magueijo & Gorski 1998; Kamionkowski & 
Jaffe 1998; Pando, Valls-Gabaud & Fan 1998; Bromley & Tegmark 
1999). An alternative approach is to analyze the present-day statis- 
tics of density or velocity fields of large scale structure (LSS). In 
principle, this is a more complicated approach, since gravitational 
instability and bias can introduce non-gaussian features in an ini- 
tially gaussian field, and these may mask the signal we desire to 
measure. Since the CMB and LSS data will shortly improve dra- 
matically with the Microwave Anisotropy Probe (MAP) and Planck 
Surveyor satellites, and the Anglo-Australian 2-Degree Field (2dF) 
and Sloan Digital Sky Survey (SDSS), respectively, it is timely to 



ask which of the CMB or LSS is the better place to look to detect a 
primordial non-gaussian signal. 

In this paper, we use the skewness and bispectrum to deter- 
mine which of the CMB or LSS provides a better probe of non- 
gaussianity. To do so, we consider several models with a primor- 
dial non-gaussianity whose amplitude can be dialed from zero (the 
gaussian limit). We then calculate the smallest non-gaussian am- 
plitude that can be detected with a plausible CMB map and with a 
plausible galaxy survey. In each case, we find that the smallest non- 
gaussian amplitude detectable with a CMB map is smaller than that 
from a galaxy survey, even if we neglect the complicating effects 
of biasing. 

Of course, there is an infinitude of possible deviations from 
gaussianity and we cannot address them all. However, physi- 
cal mechanisms that produce non-gaussianity generically produce 
a non-vanishing bispectrum. Since this is the lowest-order non- 
gaussian statistic, it is usually the most easily detectable. Although 
this argument is not fully general, we show that it holds in several 
non-gaussian models that have been considered in the recent liter- 
ature. Our results suggest that if the CMB maps provided by MAP 
and Planck are consistent with gaussian initial conditions, then any 
signatures of non-gaussianity found in galaxy surveys (apart from 
those from non-linear clustering) can be attributed to biasing. 



2 THE SKEWNESS 

In order to determine whether a field is gaussian, we need a discrim- 
inating statistic. We shall principally be concerned with the bispec- 
trum, as it is the lowest-order non-gaussian statistic that generically 
arises in physical mechanisms that produce non-gaussianity, and it 
is able in principle to distinguish between various sources of non- 
gaussianity (e.g. primordial, non-linear growth, bias). However, we 
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begin by discussing the skewness in LSS and the CMB, as this is 
simpler, and illustrates some of the effects. 

2.1 Skewness in Large scale structure 

The statistical properties of the fluctuations in the cosmological 
mass density field (5(x) = [/3(x) — p]/p can be characterized by 
the n-point moments, {(5"). By definition, (5(x)) = 0. If the fluc- 
tuation field is gaussian, then the probability distribution for 5 is 



pW = 



27ro- 



exp 



2cr2 



(1) 



from which the moments (n = 0, 1, . . .) can be calculated to be 

{<5'") = (2n - 1)!!{5')" = (2n - l)!!a'", (2) 

where a^ = {5^)- The odd moments are of course zero. To linear 
order in perturbation theory, 5 grows by an overall normalization 
factor, (5(x, t) = D{t)5{Ti, to), so an initially gaussian distribution 
will remain gaussian as long as linear perturbation theory holds. 

To higher order in perturbation theory, gravitational instabil- 
ity will induce departures from gaussianity. To describe the evolu- 
tion of non-linear fluctuations in perturbation theory we expand the 
fluctuation field in a series. 



5(i'+5(^)+. 



(3) 



where 5("'~0((5") (e.g. Goroff et al. 1986; Fry & Scherrer 1994). 
In the weakly non-linear regime, the series can be truncated to sec- 
ond order. The lowest-order deviation from gaussianity is described 
by {5'^). In the weakly non-linear regime, gaussian initial condi- 
tions give rise to a non-vanishing skewness (e.g. Peebles 1980), 



S3 



(52)2 



34 

y 



(4) 



in second-order perturbation theory (20PT). 

For generic non-gaussian initial conditions. Fry & Scherrer 
(1994) found the skewness in second-order perturbation theory 
(20PT) to be 

Q^c; , 34 26 {S'i)^ 8 3, /K(3)o] 10 C(4)o 6 /K(4)o] 
03 03,0+^ 21(5(2)3 7^°'lMW^{Si}^ 7 {S^y^' 

where the subscript denotes the quantity linearly evolved from 
the initial density field [e.g., in an Einstein-de Sitter Universe, So = 
5(zi)(l+Zi)] to the present epoch. The quantity J[^(„)o] (n — 3,4) 
is an integral that depends on the specific linearly-evolved con- 
nected (or irreducible) three- and four-point functions, C(3)o snd 
5(4)0, respectively. (For gaussian initial conditions, ^(„)o = for 
n > 3.) Fry & Scherrer (1994) find that /[C(„)o] < C(7i)o(0) for 
n = 3, 4 for several non-gaussian models they explore. All terms 
in the RHS of equation (H) are time-independent apart from ^3,0 
which scales like 53,0(2) oc £'3,0(2 = 0)(1 + z) in an Einstein-de 
Sitter Universe. 

It is therefore useful to define the time-independent quantities 
P3 = i(3),o/'^''^ (normalized skewness) and p4 = C(4),o/''"* (nor- 
malized kurtosis). When written it terms of the relevant quantities 
at decoupling (i.e., at 2 ~ 1100), 

34 



S3 



+ 



P3 



7 1100 a(z= 1100) 



+ dlP3 + d2P4, 



(6) 



where di ~ d2 — 2. 

Now suppose we have a survey of A^ independent volumes in 
which the rms fractional density contrast is a. Then, from equa- 
tions (0) and M), it follows that in the mildly non-linear regime, the 



standard error due to cosmic variance with which the skewness can 
be recovered is 



A^ 



15 

^2 



I7S32, 



(V) 



where we have used A{S^} = i/lSa^ + 10(53)2 and A(5^) ~ 
V2a^. Of course, since the mass will be traced by discrete objects 
(i.e., galaxies), the shot noise may increase the error estimate in 
equation uA. 

So now let us consider the 2dF and/or SDSS. The present-day 
skewness in volumes of side 10 h~^ Mpc could be measured with 
a standard error at least AS3 ~ 20 iV"^''^ ~ 10"\ where N is 
the number of such cubes in the survey volume. Thus, equation (p) 
tells us that a primordial normalized skewness on the 10 h^^ Mpc 
scale could be identified in a statistically-significant manner in 2dF 
and/or SDSS only if p3 exceeded about 10~^. 



2.2 Skewness in the CMB 

For simplicity, consider an Einstein-de Sitter Universe. Then a re- 
gion of comoving size lOh^^ Mpc (where h is the Hubble param- 
eter in units of 100 km sec^^ Mpc~^) subtends an angle of 0.1°. 
Now suppose that a full-sky cosmic-variance-limited 0.1 -degree- 
resolution CMB map (close to the Planck Surveyor's specifica- 
tions) finds that the distribution of temperature fluctuations AT/T 



is consistent with gaussian with a variance a 



{{Sp/pY) 



0.1{ST/Ty 



(10- 



(where 5p/p is the fractional density 



perturbation at a redshift z ~ 1100, when the CMB decouples). 
The largest primordial normalized skewness ps = {{5 p / pf) / a''' 
that would be consistent with such a map would be yJlK/Npi^ ^ 
10^"^, where A'^pix ~ 10^ is the number of 0.1 x 0.1-degree pixels. 
We have neglected instrumental noise and assumed systematic ef- 
fects will be under control. Still, a sensitivity to a value as small as 
P3 ^ 10^3, or perhaps an order of magnitude larger, is a realistic 
expectation of CMB maps. 

If we compare this with the nominal smallest normalized 
skewness ps ~ 10~^ accessible with LSS, it appears that the CMB 
has, perhaps an extra order of magnitude in sensitivity. On the other 
hand, there may be systematic effects in both the CMB and LSS 
measurements that may affect both estimates, and our argument 
was only qualitative. Thus, we conclude from this exercise that the 
CMB and LSS should provide roughly comparable sensitivity to a 
primordial skewness on 10 h~^ Mpc scales, with perhaps a slight 
edge to the CMB. Since these heuristic arguments are inconclusive 
to orders of magnitude as to which of the CMB and LSS provides 
a better probe of primordial non-gaussianity, we now proceed to 
consider realistic models more carefully. 



3 CMB AND LSS BISPECTRA 

The skewness has the advantage of being far easier to calculate 
than the full three-point correlation function, but it does not con- 
tain as much information. We therefore prefer to investigate the 
bispectrum (the three-point function in Fourier space) for reasons 
which have been rehearsed before (e.g., Matarrese, Verde & Heav- 
ens 1997). In addition, various theoretical models for structure for- 
mation yield directly Fourier-space quantities, so the bispectrum 
allows a more straightforward relation between measurable quanti- 
ties and theoretical predictions. 
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3.1 LSS Bispectrum 

We define the Fourier transform of the fractional overdensity per- 
turbation by (5k = / li^i 5(x) exp(— ik ■ x). The spatial bispec- 
trum _B(ki, k2, kg) is defined by 



{<5ki<5k,<5k3) = B(ki,k2,k3)<5°(ki +k2 +k3 



(8) 



where the angle brackets denotes an ensemble average or, by the 
ergodic theorem, the average over a large volume (fair sample), and 
(5^ is the Dirac delta function. 

To second-order in perturbation theory, the bispectrum may be 
written as a sum of a primordial part and a part induced by gravita- 
tional instability: 

B(ki,k2,k3)~Bo(ki,k2,k3) 
+ [2J(ki,k2)Po(fci)Po(fo) + cyc.] 
+ /d^k„[J(k„,k3-ka)T'=(ka,k3-ka,ki,k2) + eye], (9) 

where J(ka,kt,) is a function almost independent of the 
non-relativistic-matter density Q,q and cosmological constant A 
(Bouchet et al. 1992; Bouchet, Hivon & Juszkiewicz 1995; 
Bemardeau 1994; Scoccimarro et al. 1998; Kamionkowski & 
Buchalter 1998). Its detailed form need not concern us here. Here, 
Bit is the primordial bispectrum (which we wish to probe) linearly 
evolved to redshift z, Pq is the power spectrum and T'^ denotes the 
connected trispectrum, which is the Fourier transform of the con- 
nected 4-point correlation function. The last two terms arise from 
non-linear gravitational instability; this depends (see e.g. Catelan 
& Moscardini 1994) on the 4-point function, which has a dis- 
connected part (present in gaussian fields), and a connected part, 
T'^ (which is zero for gaussian fields). In principle, this last term 
(= Bt) may be important, as it grows as fast as the usual discon- 
nected part. We will show later that this term is very small for a 
range of proposed models. 

As already seen for the skewness, in an Einstein-de Sitter Uni- 
verse, Bo{z) ex (1 + z)~^ and [Po{z)f ex (1 + 2)"*, so the pri- 
mordial bispectrum redshifts away in comparison with that arising 
from non-linear gravitational clustering n. 

If we suppose that the galaxies (subscript g) are locally bi- 
ased with respect to the mass, then the density of galaxies Sg can 
be Taylor expanded in the mass density to second order, Sg = 
bo + biS + ^b2S , for some biasing coefficients bt. The quantity 
&o affects only k = 0, so it can be ignored. The bispectrum for the 
galaxies is then 

B9(ki,k2,k3)~6?Bo(ki,k2,k3)+ 
{Poiki)Poik2) [6?2J(ki,k2) + 626?]}+cyc. + 6?BT. (10) 

So, assuming that the initial (or the linear) power spectrum is pre- 
cisely known, we can do a likelihood analysis only if we also have 
a model for the initial (or linear) bispectrum and trispectrum as a 
function of the three wave numbers. As we found previously, the 
primordial bispectrum redshifts away, and we can realistically only 
detect non-gaussianity if Bt is significant. In section (4.1) we will 
quantify further this argument by showing that for all the models 
considered Bt <^ Bq. Let us parameterize the observed bispec- 
trum as 



* In a non-Einstein-de-Sitter model, the suppression factor (1 + 2) is re- 
placed by g(z){l + z), but for reasonable values of Qq and A, this factor 
is 0.3 < g{z) < 1. Thus, we will assume an Einstein-de Sitter model from 
now, noting that the general arguments are essentially unchanged in more 
general models. 



B5(ki,k2,k2) = Po(fci)Po(fc2)[ci2J(ki,k2) + C2]+ eye. (11) 

In order to asses the precision with which the bias could be mea- 
sured, Matarrese, Verde & Heavens (1997) evaluated the precision 
with which the parameters ei and e2 could be recovered from a 
likelihood analysis of the 2dF and/or SDSS. From their Figure 7, 
we conclude that if the positions of all of the galaxies in the survey 
volume are known, then (a) ei can be determined with an error of 
~ 6 X 10^'^ if e2 is fixed; (b) e2 can be determined with an error 
of ~ 1 X 10^^ if ci is fixed; and (c) the joint determination of 
the two parameters allows ei and C2 to be recovered with errors of 
1 X 10^^ and 4 x 10^^, respectively. The analysis of Matarrese, 
Verde & Heavens (1997) further shows that the bispectrum signal 
comes primarily from fc ~ 0.1 — 1 /i Mpc^^. (All of these esti- 
mates would be reduced only by about a factor 20 if we could map 
the mass throughout the entire Hubble volume rather than just in 
the survey volume.) We will use these results to asses the smallest 
detectable primordial non-gaussianity. 



3.2 CMB Bispectrum 

A CMB map of the temperature r(n) as a function of position fi 
on the sky can be decomposed into spherical harmonics. 



Ar(n) 
T 



y ^aimYim{n), 



(12) 



where the multipole coefficients are given by the inverse transfor- 
mation. 



a,,„ = / dnYi„,(n) — - — . 

The CMB bispectrum Bi-^i^i^ is defined by 



/ m-i mo m^\ t-> / 



^1 h I's 

mi 7712 7713 



(13) 



(14) 



where the factorization ensures statistical isotropy, and the last term 
is the Wigner 3 J symbol. 

If we parametrize the bispectrum by X times some fixed 
function — as in equation (14) — of li,l2,h, then the error on X 
is given by (see, e.g., Matarrese, Verde Sl Heavens 1997) 






E, 



y 



h I2 h 
77ii 7772 rns 

where C denotes the likelihood function and C; = (| a[" |), is the 
power spectrum of the sky fluctuation. Here we assumed that the 
departures from gaussianity are small and therefore the covariance 
matrix can be approximated by the covariance of a gaussian field 
that has the same power spectrum (Jungman et al. 1996; Heavens 
1998), and we ignore the mixing which arises from partial sky cov- 
erage. 



If one considers only the real part of {a^ 



mi m2 m^ 



) (see 



Matarrese, Verde & Heavens 1997) then n = 1/2. The quan- 
tity N[mi,li) is related to the number of non-zero terms like 
CjXJi^Ci^ in the covariance and ranges from 1 to 30. Equation 
(jlSp is valid as long as the noise does not dominate the signal. If 
A*' = 1, then the sum over the Wigner 3J symbols is unity; the fact 
that A'^ is not equal to one reduces the sum by a only few percent, 
so 
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^x 



ii<t.,<t- 



i.Bi,i 



Ci^Ci^Ci^ 



(16) 



1 



This expression ignores tlie effects of pixel noise. However, the re- 
sults we show below use I < 100 only, and for these modes, the 
pixel-noise variance in MAP and Planck will be negligible com- 
pared to the cosmic variance. 



4 SOME NON-GAUSSIAN MODELS 

We now proceed to consider several classes of physically-motivated 
models with primordial non-gaussianity in order to investigate 
more precisely the relative sensitivities of the CMB and LSS. Since 
current measurements of CMB and LSS power spectra are roughly 
compatible with cold-dark-matter (CDM) models for structure for- 
mation, we first consider CDM-like models and introduce some 
non-gaussianity in several different ways. Specifically, we consider 
models in which the gravitational potential contains a part that is 
the square of a gaussian random field and models in which the 
density contains a part that is the square of a gaussian random 
field. We note that such non-gaussianity may arise in slow-roll 
and/or nonstandard (e.g., two-field) inflation models (Luo 1994; 
Falk, Rangarajan & Srednicki 1993; Gangui et al. 1994; Fan & 
Bardeen 1992). Moreover both of these models may be consid- 
ered as Taylor expansions of more general fields, and are thus a 
fairly generic form of non-gaussianity. We also consider 0{N)-a 
models, as these will approximate the non-gaussianity expected in 
topological-defect models. 

4.1 Quadratic model for the potential 

We start by considering a model in which the gravitational potential 
$ (in conformal Newtonian gauge) is a linear combination of a 
gaussian random field (j> and a term proportional to the square of 
the same random field. 



<i> = cb + a{cb^-{4>^}), 



(17) 



where a parametrizes the non-gaussianity; in the limit a ^ 0, 
the model becomes gaussian.P] These models contain a primordial 
bispectrum for the gravitational potential. 



B*(ki,k2,k3) ~ 2Q[P<l,(fci)P<l,(fc2) +CJ/C.]. 



(18) 



The leading terms in the connected trispectrum for the $ field will 
be 



T|(ki 



■,k4) 



4a2Ps(ki)P<s,(k2)[P*(lki + k3|)+P*(|ki+k4l)] + cyc. (19) 

We use a scale-invariant primordial spectrum, Pj. = Ank'^^, 
where the amplitude is fixed by COBE to be, in our Fourier trans- 
form conventions. Ah — 10^^". 

The relative contribution to the bispectrum of Bq and Bt, 
R = Bq/Bt, can be obtained from equation (9) and (19). By eval- 
uating the integral we find: 

^' (20) 

oAh 



R 



T We note that such non-gaussianity can arise in standard slow-roll infla- 
tion, and the parameter a can be related to inflaton-potential parameters 
(e.g., Falk, Rangarajan & Srednicki 1993; Gangui et al. 1994; Wang & 
Kamionkowski 1999). 
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Figure 1. The inverse of the function A^ j. evaluated at redshift z = 0. In 
the regime where 20PT holds, this coefficient is between 10^^ and 10^^. 



Thus, Bt ^ Bo for Q ^ 10^", so we can safely neglect the contri- 
bution from the trispectrum. The Bt contribution to the bispectrum 
will therefore be negligible compared with Bo if a <C {Ah)~^- 

The Fourier coefficients <l?(k) for the gravitational poten- 
tial remain constant to linear order in perturbation theory in an 
Einstein-de- Sitter Universe. The Fourier coefficients S(k, z) of the 
density field evolve with time and are related to those of the gravi- 
tational potential by the Poisson equation, which can be written 



S{k,z) = Mkiz)${k), wheieMkiz) 



2k^T{k){l + z) 



3Hi 



(21) 



The inverse of the function A4k{z = 0) is plotted as a function of 
wavenumber k in Fig. Ill Thus, the linearly-evolved power spectrum 
for the mass, Po(fc, z) at redshift z is related to that for the gravi- 
tational potential (which remains constant in and Einstein-de-Sitter 
Universe) by 

Po{k,z)^[Mk{z)fP^{k), (22) 

where the transfer function is (e.g., Bardeen et al. 1986) 

T{k) ^ ,^ , „,_ , i,,,,, , ^,_,, , (23) 



(l + Bfc + Cfc»/2 +Dfc2)' 

and B = 1.7Mpc(J7^^)-\ C = W\ife'l'^{Q.h^Y'^^'^ and D = 
IMpc^ {QJn}')~'^ . The bispectrum for the mass for this model is thus 



B(ki^ 

Ml 



k2,k3 



{p(fc.)P(fe2) [(2a^,^) 



+ 2J(ki,k2 



+ eye. (24) 



Comparing equation ( |24| ) with equation (Ull), we see that this 
particular form of primordial non-gaussianity leads to a present- 
day bispectrum that looks like a scale-dependent non-linear bias. 
In particular, for equilateral configurations, equation (E4) becomes 
identical to equation mm if we identify C2 = IolJ M^. We can 
therefore use the results of Matarrese, Verde & Heavens (1997) to 
determine the smallest a that would be detectable by 2dF/SDSS, 
under the assumption that there were no non-linear biasing. In the 
range where most of the signal comes from, M^^ — 10~^ — 10~^. 
Thus, we conclude that the smallest a that would give rise to an 
observable signal in the 2dF/SDSS bispectrum is a ~ 10^ — lO**. 

The primordial gravitational-potential bispectrum in equation 
( |l8| ) will lead to a non-zero bispectrum in the CMB via the Sachs- 
Wolfe effect, and this bispectrum can be calculated to be (Luo 



© 0000 RAS, MNRAS 000, 000-000 



Large-scale structure, the cosmic microwave background, and primordial non-gaussianity 



1994; Wang & Kamionkowski 1999), 



t 



(2;i + l)(2/2 + l)(2/3 + l) 
47r 
2a 



h 




h 




j4sw 



[Ci-^^Ci2 + Ci-^Cig + Ci^Ci.:^] , (25) 



where Asw — 1/3 is the Sach-Wolfe coefficient. Plugging this 
bispectrum (and a scale-invariant set of Ci) into equation (Ua), we 
learn that the smallest a that could be detectable with a CMB map 
(using only I < 100) is ~ 20fI (The dashed curve in Fig. shows 
the smallest a detectable with the CMB as a function of the largest 
multipole moment I used in the analysis.) Thus, we conclude that 
the CMB will be at least two orders of magnitude more sensitive to 
a non-zero value of a than LSS. 




4.2 Quadratic model for the density 

We now consider an alternative model in which the density field 
(rather than the gravitational potential) contains a term that is the 
square of a gaussian random field. 



5 = 



<P + a {<P'' - {^'')) . 



(26) 



where now is some other gaussian random field. Such a model 
has been considered in some two-field inflation models (e.g. Luo & 
Schramm 1993). 

Since the density perturbation evolves in linear theory with 
time (unlike the gravitational potential), we must specify the epoch 
at which the density perturbation S is related to (j> through equa- 
tion (pfil). We choose this epoch to be the current epoch, z = 0. 
Doing so, the spatial mass bispectrum today due to primordial non- 
gaussianity is 



Bo(ki,k2,k3) ~ 2aP{ki)P{k2) + eye, 



(27) 



and Luo & Schramm (1993) have shown that the contribution to 
the present-day bispectrum from the primordial trispectrum is neg- 
ligible compared with this (Bt <^ Bo), as long as 20PT holds. By 
comparing with equation (|ll[), we see that this form of primordial 
non-gaussianity gives rises to a present-day skewness that mimics 
precisely that due to a scale-independent non-linear bias. Again, 
from the results of Matarrese, Verde & Heavens (1997), the small- 
est detectable a in this model is ^ O.Ol. 

Now let us consider the CMB bispectrum of this model. The 
spatial bispectrum for the gravitational potential here is 

2aMkiMk2 



B$(ki,k2,k3 



M, 



[P<i.(fci)P*(fc2)] + cyc. (28) 



Although we have not carried out an exact calculation of the bis- 
pectrum for this model, it is easily seen (at least in the / 3> 1 limit) 
to be 



(2;i + l)(2/2 + l)(2/3 + l) 



47r 



h 




^3 





■f Actually, our calculation ignores the physics that gives rise to the acoustic 
peaks at I > 100. However, we restrict our analysis to I < 100. A more 
accurate treatment would change our results by no more than an order of 
magnitude, and this level of precision is sufficient for our purposes. 
3 This result is in rough agreement with the conclusion that the o ^ 10 — 
100 (Luo 1994; Gangui et al. 1994; Mollerach et al. 1995) non-gaussian 
signal from non-linear evolution prior to ^ ~ 1100 is at best marginally 
detectable (Spergel & Goldberg 1999). 



Figure 2. The smallest non-gaussian amplitude a that can be detected with 
the CMB as a function of the largest multipole moment I used in the anal- 
ysis. The solid line is for a model in which the present-day fractional over- 
density field contains a term that is the squai'e of a gaussian random field; 
the dashed line is for a model in which the gravitational potential contains 
a term that is the square of gaussian random field; and the dotted line refers 
to the 0{N) a models. 
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Applying equation (llq), we find that the smallest detectable a in 
this model is ~ 0.01 (using only I < 100; results for other Zmax are 
shown in Fig. bl), which is comparable to the LSS error. However, 
since this model has more power in the CMB bispectrum at larger 
I, the smallest detectable a decreases by roughly an order of mag- 
nitude even if we go out only to /max ~ 200. Also in this case we 
find that CMB will provide a more precise probe of a primordial 
bispectrum. Moreover, as a corollary, this particular result demon- 
strates that if the maps provided by MAP and Planck are consistent 
with gaussian, then any measurement of a non-zero C2 from the 
LSS bispectrum can be interpreted unambiguously as evidence for 
non-linear biasing, rather than some primordial non-gaussianity. 



4.3 0{N) a models 

The 0{N) sigma model provides an approximation to the non- 
gaussianity expected in topological-defect or scalar-field-alignment 
models (Turok & Spergel 1991; Jaffe (1994)). The TV = 1 
model has domain walls, N — 2 global strings, A'^ = 3 global 
monopoles, A^ = 4 global textures, and higher A'^ correspond to 
non-topological-defect models. For large A'^ this model approaches 
the gaussian model, so we take a — N~^'^ and as a — > 
the models become asymptotically gaussian. Since calculation of 
power spectra and higher-order statistics for the CMB and LSS are 
quite involved for these models, our analysis will be only approx- 
imate. As we will see below, these order-of-magnitude estimates 
will be sufficiently precise for our purposes. 

For equilateral-triangle configurations and for values of k that 
can be probed with 2dF/SDSS, the power spectrum, bispectrum and 
trispectrum for LSS in the linear regime are (Jaffe 1994) 



P~12.5K:'^kT\k), Bo~l.6K.^aT^{k), 






(30) 



(31) 
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where /C ~ 30- 100 Mpc^ h'^. From equations @ and ^, Bt 
is found to be always < 0.4i3o- Since Bq will be the dominant con- 
tribution to the bispectrum we conclude that the non-gaussianity of 
this model will act like a scale-dependent non-linear-bias contribu- 
tion with 

C2 ~ ^ ~ a[ICk^lO''T{k)]-\ (32) 

For the scales probed by SDSS/2dF, C2 < q/300, and therefore the 
minimum a detectable from LSS will be ~ 30. 

Precise calculation of the CMB bispectrum for these mod- 
els is well beyond the scope of this paper To obtain an order- 
of-magnitude estimate, we assume that the spatial polyspectra for 
this model, equations (pO|)-(pl[), give rise to potential polyspectra 
through the Poisson equation and then that the CMB anisotropy is 
proportional to the potential perturbation at the surface of last scat- 
ter. Doing so, we finqj 



{2li + I){2l2 + I){2l3 + 1) 



An 



h 




^3 





x^QQQ a[CnCi^Ci^f''-^ (33) 



Using equation (|16[), we find that the smallest detectable a would 
be ~ 0.1 (using multipole moments up to /max = 100; the depen- 
dence on Zmax is indicated in Fig. ha). Thus, we again conclude 
that the CMB will provide a more precise probe of a primordial 
non-gaussianity. 



5 CONCLUSIONS 

We addressed the question of which of the CMB or LSS is bet- 
ter poised to detect primordial non-gaussianity of several varieties. 
We used the bispectrum as a discriminating statistic since it is the 
lowest-order quantity that has zero expectation value for a gaussian 
field. We considered three forms of non-gaussianity: one in which 
the gravitational potential contained a term that was the square of a 
gaussian field; another in which the density field was the square of 
a gaussian field; and a third that resembles that expected from topo- 
logical defects. We showed that in all cases, the CMB is likely to 
provide a better probe of such non-gaussianity. One of these models 
produced a mass bispectrum that mimicked a scale-dependent non- 
linear bias, and the others mimicked a scale-independent non-linear 
bias. 

Of course, our results are not fully general. In principle, it is 
possible to think of some other type of non-gaussianity for which 
our conclusions would not hold. However, plausible physical mech- 
anisms that produce nearly scale-invariant power spectra should 
generically produce non-gaussian signals that have scale depen- 
dences roughly like those that we investigated. Thus we may con- 
clude that if CMB maps turn out to be consistent with gaussian 
initial conditions, any non-gaussianity seen in the LSS bispectrum 
can be unambiguously attributed to the effects of biasing. 
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